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RAPID COOLING OF DUSTY GAS IN ELLIPTICAL GALAXIES 
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ABSTRACT 

We propose a stellar origin for the central dust clouds observed in most giant elliptical galaxies. Dusty 
gas ejected from evolving red giant stars in E or cD galaxies can cool rapidly even after entering the 
hot, X-ray emitting gas. Cooling by thermal collisions with dust grains can be faster than either the 
dynamical time in the galactic potential or the grain sputtering time. Some grains survive in the cooled 
gas. Dusty stellar outflows cool more efficiently in the central regions where the stellar metallicity is 
higher. Mergers with gas and dust-rich dwarf galaxies may occasionally occur but are not required to 
explain the observed dust clouds. 

Subject headings: galaxies: elliptical and lenticular, CD - ISM: dust - cooling flows - X-rays: galaxies - 
galaxies: clusters: general - X-rays: galaxies: clusters 



1. INTRODUCTION 

While most of the interstellar gas mass in elliptical 
galaxies, typically ~ 10 9 Mq, is at the galactic virial tem- 
perature T V i r ~ I0 7 K, the cores of ~ 80 percent of all 
large E galaxies contain dusty clouds of cold gas of mass 
< I0 4 - 10 5 M Q (e.g. van Dokkum & Franx 1995). In 
addition, diffuse optical line emission from gas at T ~ 10 
K is observed in most or all elliptical galaxies within a few 
hundred parsecs from the center (e.g. Caon et al. 2000), 
particularly in those galaxies formerly thought to contain 
"cooling flows." The notion that the hot gas cools to low 
temperatures is not supported by XMM spectra (e.g. Pe- 
terson et al. 2001; Xu et al. 2001), i.e., the observed 
cooling rate M is lower than that expected in traditional 
cooling flows, ~ 1 Mq yr _1 . In any case, cooling of es- 
sentially dust-free hot gas cannot produce dust clouds. In 
the following we propose that the dusty clouds are formed 
from gas expelled from evolving metal-rich red giant stars 
in E galaxy cores. Even if these dusty stellar envelopes are 
heated to ~ T V i r after leaving the stars, the gas can cool 
rapidly by collisions between thermal electrons and dust 
grains. 

2. COOLING-SPUTTERING EVOLUTION 

Direct evidence that mass is being lost from evolving 
stars in X-ray luminous ellipticals is provided by near-IR 
(~ 10/xm) observations of hot circumstellar dust (Knapp 
et al. 1989, 1992; Athey et al. 2002). The 9.7 fiia peak 
from (oxygen-rich) silicate grains apparent in ISO obser- 
vations places an upper limit on the grain radius, a < 1 
(im (Figure 5 of Laor & Draine 1993). According to Athey 
et al., the total luminosity of excess near-IR radiation is 
consistent with expected total stellar mass loss rates in 
giant E galaxies, ~ 1 Mq yr _1 . 

However, as dusty gas moves away from orbiting, mass- 
losing red giant stars, it encounters the hot interstellar 
gas at T ~ T v i r ~ 1 keV and is violently disrupted by 
hydrodynamic instabilities. The interaction of the stellar 
wind/planetary nebula with the hot gas is extremely com- 



plicated, but this ejected gas cannot remain at T ~ 10 4 
K for more than tn a ~ 10 5 years without exceeding 
the typical Ka + [Nil] luminosities observed (Mathews 
& Brighenti 1999). XMM spectra show no evidence for 
slowly heating gas at T < T V i r /2. This empirical evi- 
dence supports rapid thermal mixing of the stellar ejecta, 
which is probably a consequence of gasdynamical insta- 
bilities that enormously increase the surface area between 
the stellar ejecta and the ambient hot gas. In view of the 
complexity of this interaction, the thermal evolution T(t) 
of gas ejected from stars is difficult to estimate, but limits 
can be set. In the following discussion we adopt a use- 
ful limiting assumption for T(t): the dusty gas is rapidly 
heated to the temperature of the local hot gas with the 
dust intact. This assumption maximizes the likelihood 
that the stellar gas will remain thermally merged with the 
hot gas phase. Grain sputtering during the rapid heating 
process is unlikely since the heating time must be ;$ tfjai 
which is very much less than the sputtering time (see be- 
low). 

When dust grains are introduced into a hot plasma with 
T > 3 x 10 6 K, they are sputtered (eroded) by impacting 
thermal H + and He +2 ions. The grain radius a decreases 
at a rate 
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cm s 
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(Draine & Salpeter 1979; Tsai & Mathews 1995) where 
n{He ++ ) = n{H+)/lQ, Here h w = 3.2 x 10~ 18 cm 4 s" 1 
and T s = 2 x 10 6 K are fitting parameters, n p — n e (4 — 
3^t) / (2+/i) is the proton density and fi = 0.61 is the molec- 
ular weight. We assume for simplicity that the cooling 
occurs at constant gas pressure P so n e = (2 + \i)Pjh\ikT. 

As soon as the dusty gas is heated (by assumption) to 
the local gas temperature, it will begin to cool. Cool- 
ing occurs by thermal X-ray emission from the gas and 
by inelastic impacts of thermal electrons with the grains. 
Grain-electron cooling can greatly exceed that from ther- 
mal X-ray emission (e.g. Dwek 1987; Dwek & Arendt 
1992) provided the grains are not sputtered away. If the 
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cooling time is shorter than the dynamical time in the 
galactic potential, the gas will tend to cool approximately 
in situ, the case we consider here. During isobaric cooling, 
the gas temperature decreases according to 



dT 
~dt 



2 TV 
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(A + A d ) 



2 P^. (2) 



Here A(T, z) is the usual thermal radiation coefficient for 
gas with metallicity z in solar units (Sutherland & Dopita 
1993) and A^ is the coefficient for dust-induced cooling. 

The cooling of the hot plasma by grains is due mostly to 
inelastic impacts of thermal electrons with the grains. For 
high temperature plasmas, such as we consider, the rate of 
electron-grain collisions is unaffected by the grain charge 
(Draine & Salpeter 1979) which we assume is regulated by 
emission of low energy electrons from the grain surface. 
An individual grain cools the local gas at a rate 

e g = A-Ka 2 {l/A)n e v e E e = 1.27 x 10~ 7 n e T 7 3/2 a 



erg s 



(3) 

where E e = (3/2)fcT and v e = (SkT/irnie) 1 / 2 are the 
mean energy and velocity of electrons, T7 = T/(10 7 K) 
and a M is the grain radius in microns. At pressures of in- 
terest, -P-io = P/(1CT 10 dy), the plasma mean free path, 
A « 0.02T 3 /P_io cm, is small; only electrons in the cool- 
ing regions collide with the grains; electrons from the hot 
ambient gas thermalize in the cooling plasma before col- 
liding with a grain. In NGC 4472, a giant elliptical in the 
Virgo cluster which we use for reference, n e « O.OSSlr^; 18 
cm" 3 so e g = 7.7 x 10~ 9 a 2 and 7.3 x lCT 10 ^ at r = 1 kpc 
and the effective radius r = R e = 8.57 kpc respectively. 
We adopt a distance of 17 Mpc to NGC 4472. 

Grains can internally absorb the energy of thermal elec- 
trons provided E e < = 23a 2 / 3 keV (Dwek & Werner 
1981), which we assume is generally satisfied. Cooling 
by electron impact on grains during the sputtering time, 
t sp w a/\da/dt\ w 1.2 x 10 8 a M (n e /10 _2 )~ 1 yrs, greatly ex- 
ceeds the additional gas cooling required to fully ionize the 
grain atoms ejected in the sputtering process or to bring 
these sputtered ions up to local thermal velocities. 

At any moment the space density of interstellar grains 
n g in luminous elliptical galaxies is very inhomogeneous, 
being concentrated near mass-losing stars. If m e j is the 
typical total mass of the stellar envelope ejected by evolv- 
ing stars and p g = 3.3 gm cm~ 3 is the density of (silicate) 
grains, then the total number of grains expelled by a single 
star is Af g — 5m e j / p g (4:/3)ira 3 where 5 is the ratio of the 
dust mass to the gas mass when ejected from the stars. We 
assume that 5 = O.Olz^Fe decreases with stellar metallic- 
ity as in NGC 4472 where z*, Fe ~ 0.675(r/ J R e )-°- 207 in 
solar units (Brighenti & Mathews 1999). 

The initial space density of grains is proportional 
to the local gas density n g0 = (Af g /m e j)p. As the 
gas cools isobarically from some initial temperature Tq, 



the grain space density increases, n g 
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[Sm p / p g (4/3)ira 3 ](pP/k)/T. The dust cooling coefficient 
in Equation 2, A^ = n g e g /(p/m p ) 2 , is then 
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To explore the influence of grain cooling we integrate 
Equations 1 and 2 beginning with gas at temperature 



To = T(r) and pressure P(r) at galactic radius r in 
the hot gaseous atmosphere of NGC 4472. The pressure 
P = n e (r)kT(r)5/(2 + p) is constant during cooling. For 
simplicity, we assume the gas is filled with grains all hav- 
ing the same initial radius do- As the gas cools the sput- 
tered debris from the grains increases the metallicity of the 
gas phase. To evaluate the plasma cooling rate coefficient 
A(T, z) during the cooling, we assume that the gas phase 
metal abundance in solar units increases with decreasing 
grain mass (oc a 3 ), z/zq = 1 — (a/ao) 3 (<5/<5 maa; ), where 
Smax ~ 0.053 is the maximum ratio of dust to gas mass if 
all elements heaver than He were in grains. 

3. RESULTS 

We describe solutions of Equations 1 and 2 for the el- 
liptical galaxy NGC 4472 with initial conditions at radii 
r = 1 and r = R e = 8.57 kpc. According to our assump- 
tion, the dusty gas expelled from a star is rather quickly 
heated to the local temperature of the hot gas at these two 
galactic radii, T = 0.95 x 10 7 and 1.16 x 10 7 K, where the 
gas density is n e0 = 6.52 x 10~ 2 and 4.45 x 10~ 3 cm~ 3 
respectively. Figures la and lb show the combined evolu- 
tion of the gas temperature T(t) and grain radius a(t) for 
grains of initial radius do = 1 and 0.1 micron, assuming 
S = 0.01 and 5 — 0.0064 at r = 1 and r = 8.57 kpc respec- 
tively. For comparison we also plot the cooling curve T(t) 
of dust-free gas (do = Op) with the same T and n e0 . 

Dust-enhanced cooling can reduce the cooling time of 
recently ejected gas by over an order of magnitude com- 
pared to the dust-free case. For fixed S, the ratio A^/A 
increases rapidly with decreasing grain radius ao because 
of the larger total surface area of smaller grains. 

A defining characteristic of cooling flows is that the iso- 
baric radiative cooling time t coo ^ ra d = 5m p kT/2ppA w 
10 8 r^. p 2 yrs (in NGC 4472) exceeds the dynamical time 

tdy ~ [r 3 /GM(r)] 1 / 2 in the galactic potential. Dust-free, 
moderately positive density perturbations do not cool ap- 
preciably faster than unperturbed regions. If overdense 
regions remain coherent, they oscillate radially in the hot 
gas atmosphere, becoming underdense relative to the am- 
bient gas during half their cycle (e.g. Loewenstein 1989). 

The dust-rich regions we describe here are not initially 
in an overdense state, but become overdense as the cool- 
ing proceeds. The dust-rich gas cools locally if the cooling 
time is less than the dynamical time at the two selected 
galactic radii, td y = 2.1 x 10 6 at r = 1 kpc and 1.2 x 10 7 
yrs at 8.57 kpc. Figures la and lb show that t coo i < td y 
is possible at r = 1 kpc, so gas at small galactic radii can 
cool to low temperature and also preserve a small amount 
of its original dust since t coo i < t sp . The cooling curves in 
Figures la and lb are independent of m e j and are there- 
fore unchanged if the dusty envelope is disrupted or frag- 
mented during its interaction with the hot gas. We show 
one case in Figure lb in which dusty gas with a = 0.1 
pm begins cooling at To/3, before fully thermalizing with 
the hot gas. The cooling time is significantly reduced, con- 
firming our assumption that T(t = 0) = T maximizes the 
cooling time. 

As dust-rich stellar winds flow away from orbiting red 
giant stars, the winds are decelerated by interactions with 
the local hot gas, but the grains, impelled by their higher 
momentum, can move into the ambient gas. It is possible 
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in principle that the ejected dust cloud occupies and cools 
a gas mass larger than that ejected from the star, m e j. If 
the cooled mass exceeds m e j by a factor u>, then the cor- 
responding cooling evolution is found simply by replacing 
S in Equation 4 with S/lo. In Figure lc we show cooling 
curves for u> = 2 at the same two galactic radii for grains of 
initial radius oq = 0.1 micron. Both cooling times exceed 
the dynamical time at the galactic radii considered, so the 
cooling gas may have a more complex dynamical evolution 
that may inhibit the cooling. 

4. DISCUSSION AND CONCLUSIONS 

Stellar mass loss has been regarded as an important in- 
ternal source of hot gas within group or cluster-centered 
E galaxies. A ~ 13 Gyr old stellar population of total 
mass M* t typically expels ~ 1.5(M* t /10 12 M ) M Q 
yr~ 4 . In cluster-centered galaxies, such as M87 in Virgo 
and NGC 4874 in the Coma cluster, the gas tempera- 
ture rapidly decreases in the central ^15 kpc from the 
virial temperature of the cluster (3-8 keV) to the stellar 
virial temperature ~ 1 keV at r < 3 kpc (Molendi 2002; 
Vikhliin ct al. 2001). Detailed gasdynamical models of 
traditional cooling inflows in these cluster-centered galax- 
ies indicate that this steep temperature gradient can only 
be understood if gas is being cooled by thermalization of 
stellar ejecta (Brighenti & Mathews 2002a). The increase 
in the hot gas oxygen abundance toward the center of M87 
(Gastaldello & Molendi 2002) is expected if mass lost from 
the SNII-enriched stars has thermally mixed into the hot 
ISM and if the stars are more O-rich than the Virgo clus- 
ter gas. Nevertheless, we propose here that the dust com- 
ponent can cool the stellar ejecta soon after it thermally 
merges with the ambient hot gas, particularly at r < 1 
kpc. 

We now review how dust-enhanced cooling applies to 
the issues discussed in the Introduction. 

(1) What is the source of interstellar dust observed in 
the cores of ~ 80 percent of all large E galaxies? van 
Dokkum & Franx (1995) and others describe small central 
dust disks, lanes or clouds typically a few 100 pc in size. 
The small masses of dust ^ 10 — 10 5 M Q in these cores 
can easily be produced in ~ 10 8 — 10 9 yrs as dust-rich 
gas cools near the galactic center with incomplete sput- 
tering. It is difficult to accurately estimate the rate that 
dust accumulates near the center without knowledge of 
the dust size distribution, the exact thermal history T(t) 
of the stellar ejecta, or the star formation rate in the cold, 
dusty clouds. Since Type la supernova remnants occupy a 
very small fraction of the interstellar volume at any time 
(Mathews 1990), heating by Type la supernovae is unlikely 
to interfere with the cooling we describe here. Heating by 
active galactic nuclei has often been suggested to explain 
why the hot gas in E galaxies fails to cool to low tempera- 
tures (e.g., Rosner & Tucker 1989; Binney & Tabor 1995; 
Brighenti & Mathews 2002b). If such heating occurs, it 
must be gentle enough not to destroy the observed central 
dust clouds or this dust must be rapidly regenerated. 

(2) What is the source of the diffuse optical line emission 
(Ha + [Nil]) observed in most or all cooling flow galax- 
ies? The velocities of this diffuse emission are unrelated 
to stellar velocities (Caon et al. 2000). In some E galax- 
ies the optical line emission spatially correlates with X-ray 



features (e.g. Trinchieri & Goudfrooij 2002) or traces the 
perimeters but not the centers of X-ray cavities (McNa- 
mara, O'Connell & Sarazin 1996; Blanton ct al. 2001). 
While these observations suggest that some gas cools to 
~ 10 4 K from the hot phase, there is no evidence for this 
in XMM X-ray spectra. However, even if dusty gas ejected 
from metal-rich stars is heated to ~ T V i r , it can quickly 
cool back to ~ 10 4 K where the cooling may be temporar- 
ily arrested by absorption of galactic UV starlight (Bincttc 
et al. 1994). Exactly what happens next is unclear, but 
the dust can help cool the gas to much lower temperatures 
if the clouds become optically thick to UV radiation. It 
has often been suggested that dusty gas at T ~ 10 4 K de- 
rives from mergers with gas-rich dwarf galaxies (e.g. Caon 
et al. 2000; Trinchieri & Goudfrooij 2002); this can be 
verified if the 10 4 K gas is counter-rotating. Nevertheless, 
cold, dusty gas can arise naturally from stars in E galaxy 
cores. 

(3) Why is the observed cooling rate M from X-ray ob- 
servations so low? Xu et al. (2001) found from RGS XMM 
observations of elliptical galaxy NGC 4636 that X-ray lines 
expected from gas cooling near T ~ 2 x 10 6 K are unusu- 
ally weak, indicating a total cooling rate of M < 0.30 M 
yr -i j n r < 2 5 kp C (f or d — 17 Mpc). FUSE observations 
of OVI lines in NGC 4636 (emitted at T ~ 3 x 10 5 K) in- 
dicate M k, 0.17±0.02 M Q yr _1 within 1.2 kpc (Bregman 
et al. 2001). Both cooling rates are less than the ^1 — 2 
M Q yr _1 predicted in traditional cooling flow models (e.g. 
Bertin & Toniazzo 1995). 

Rapid dust-assisted cooling may be relevant to this dis- 
crepancy in two ways: (1) by reducing the total rate that 
gas enters the hot phase, and (2) by reducing the X-ray 
emission from cooling thermal gas at subvirial tempera- 
tures. In models of traditional cooling flows, gas ejected 
from stars is usually assumed to enter the hot interstellar 
gas. The stellar mass loss rate in a giant E galaxy, ~ 1M® 
yr~ 4 , is an important source of gas since it is comparable 
to the expected cooling rate of the hot gas. In gas dy- 
namical models the cooling rate M near the center of the 
flow varies inversely with the specific rate of stellar mass 
loss, a* = M*/M*. This is true even when there is an 
extended reservoir of circumgalactic hot gas due to cosmic 
accretion onto the surrounding galaxy group. Therefore, 
the apparent cooling rate M observed with XMM should 
be reduced approximately in proportion to the fraction of 
stellar gas that fails to enter the hot gas because of dust- 
assisted cooling. However, if more than <~ 90 percent of 
the stellar ejecta fails to enter the hot phase, our gasdy- 
namical models indicate that the density of hot interstel- 
lar gas is lowered sufficiently to initiate a strong galactic 
wind driven by Type la energy. The transition to wind 
flows is rather sudden. Since the very low L x characteris- 
tic of strong winds are not observed, this sets a limit on 
the efficiency of dust-assisted cooling described here. 

If dusty gas from stars is indeed heated to ~ T V i r , XMM 
observations require that it not emit thermal X-rays at 
intermediate temperatures either as it is heated or as it 
cools afterward. We argue above that the heating phase 
is likely to be rapid, ^ tu 

a ~ 10 yrs. When heated to 
~ T v i r , the stellar gas may undergo rapid dust-assisted 
cooling, during which the X-ray emission from intermedi- 
ate temperatures is greatly reduced. The radiation ex- 



4 



pected from an X-ray line of emissivity su ne (T, p) ergs 
s -1 gm _1 , E Une = J e Hne dT/\dT/dt\ ergs gm -1 , varies 
inversely with the cooling rate \dT/dt\, which can be > 10 
times larger than normal plasma cooling. The X-ray line 
emission is suppressed by a similar factor. 

In dust-assisted cooling, much of the thermal energy in 
the hot gas is radiated by dust in the far-IR. The maxi- 
mum FIR luminosity expected, Lpm ~ 5MkT/2/j,m p ~ 
2 x 10 41 (T/10 7 )(M/M Q yr -i) erg s" 1 , is consistent with 
ISO FIR luminosities observed in NGC 4472 and NGC 
4636, fL„<2-5x 10 41 erg s" 1 (Tcmi ct al. 2003). 

Nevertheless, it seems unlikely that dust-enhanced cool- 
ing can explain the weakness of X-ray cooling lines from 



gas at subvirial temperatures in clusters (e.g. Peterson 
et al. 2001) where the expected cooling rates, M > 100 
Mq yr -1 , greatly exceed the stellar mass loss rate in the 
central E or cD galaxy 
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Fig. 1. — (a): Evolution of the gas temperature (solid lines) and grain radius in fi (short dashed lines) for dust-induced cooling in NGC 
4472 at the effective radius, r = 8.57 kpc. Each set of curves is labeled with the initial grain radius ao in /im and the solid curve for ao = 
is the temperature evolution for dust-free gas at the same radius. The ratio of the rates of dust and gas cooling A^/A is shown with long 
dashed lines. The vertical heavy line shows the dynamical time tjj for freefall from r to the galactic center, (b): Same evolution as in (a) but 
at radius r = 1 kpc in NGC 4472. The heavy lines show the T, a-evolution of a = 0.1 /i grains if cooling commences at a lower temperature 
To/3, (c): Two sets of evolutionary curves when the grains are dispersed into a mass of gas equal to twice that ejected from the star. Both 
sets of curves are for ao =0.1 (im and are labeled with the initial galactic radius. 
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